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Abstract 

The main purpose of this paper is to explore the solutions of 
Bianchi type III cosmological model in Brans Dicke theory of gravity 
in the background of anisotropic dark energy. We use the assumption 
of constant deceleration parameter and power law relation between 
scalar field <f> and scale factor a to find the solutions. The physi- 
cal behavior of the solutions has been discussed using some physical 
quantities. 
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1 Introduction 

Recent experimental data [1] about late time accelerated expansion of the 
universe has attracted much attention in the recent years. Cosmic accelera- 
tion can be well explained from high red-shift supernova experiments . The 
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recent results from cosmic microwave background fluctuations [2] and large 
scale structure [Hj suggest the expansion of universe. Dark energy seems to 
be best candidate to explain cosmic acceleration. It is now believed that 96 
percent energy of the universe consist of dark energy and dark matter (76 
percent dark energy and 20 percent dark matter) [2], H] . Dark energy is the 
most popular way to explain recent observations that the universe is expand- 
ing at an accelerating rate. The exact nature of the dark energy is a matter of 
speculation. It is known to be very homogeneous, not very dense and is not 
known to interact through any of the fundamental forces other than gravity. 
Since it is not very dense, roughly 10 — 29 grams per cubic centimeter, so it is 
a difficult task to detect it in the laboratory. It is thought that dark energy 
have a strong negative pressure in order to explain the observed acceleration 
in the expansion rate of the universe. 

Dark energy models have significant importance now as far as theoretical 
study of the universe is concerned. It would be more interesting to study 
the variable equation of state (EoS), i.e. P = pco(t), where P is the pressure 
and p is the energy density of universe. Usually EoS parameter is assumed 
to be a constant with the values —1, 0, — | and +1 for vacuum, dust, 
radiation and stiff matter dominated universe, respectively. However, it is a 
function of time or redshift [5] in general. Latest observations [6] from SNe 
la data indicate that u is not constant. In recent years, many authors [7]-|ll] 
have shown keen interests in studying the universe with variable EoS. Sharif 
and Zubair [7] discussed the dynamics of Bianchi type VIq universe with 
anisotropic dark energy in the presence of electromagnetic field. The same 
authors [8] explored Bianchi type / universe in the presence of magnetized 
anisotropic dark energy with variable EoS parameter. Akarsu and Kilinc j9] 
investigated the general form of the anisotropy parameter of the expansion 
for Bianchi type III model. 

The isotropic models are considered to be most suitable to study large 
scale structure of the universe. However, it is believed that the early uni- 
verse may not have been exactly uniform. This prediction motivates us to 
describe the early stages of the universe with the models having anisotropic 
background. Thus, it would be worthwhile to explore anisotropic dark en- 
ergy models in the context of modified theories of gravity. Among the var- 
ious modifications of general relativity (GR), the Brans-Dicke (BD) theory 
of gravity [12] is a well known example of a scalar tensor theory in which the 
gravitational interaction involves a scalar field and the metric tensor. One 
extra parameter w is used in this theory which satisfies the equation given 
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by 

8ttT 

where is known as BD scalar field while T is the trace of the matter 
energy-momentum tensor. It is mentioned here that the general relativity 
is recovered in the limiting case w — > oo. Thus we can compare our results 
with experimental tests for significantly large value of w. 

Bianchi type models are among the simplest models with anisotropic 
background. Many authors [13]- [23] explored Bianchi type spacetimes in 
different contexts. Moussiaux et al. [25] investigated the exact solution for 
vacuum Bianchi type-Ill model with a cosmological constant. Xing-Xiang 
[2"E] discussed Bianchi type III string cosmology with bulk viscosity. He 
assumed that the expansion scalar is proportional to the shear scalar to find 
the solutions. Wang [27] investigated string cosmological models with bulk 
viscosity in Kantowski-Sachs spacetime. Upadhaya [28] explored some mag- 
netized Bianchi type-/// massive string cosmological models in GR. Hellaby 
[2U] gave an overview of some recent developments in inhomogeneous models 
and it was concluded that the universe is inhomogeneous on many scales. 
Sharif and Shamir [T7J, [18] have studied the solutions of Bianchi types I and 
V spacetimes in the framework of f(R) gravity. Recently, we JT9J [20] ex- 
plored the exact vacuum solutions of Bianchi type J, III and Kantowski- 
Sachs spacetimes in the metric version of f(R) gravity. 

The study of Bianchi type models in the context of BD theory has at- 
tracted many authors in the recent years [30]. A detailed discussion of BD 
cosmology is given by Singh and Rai [3T] . Lorenz-Petzold [32] studied exact 
Bianchi type-Ill solutions in the presence of electromagnetic field. Kumar 
and Singh [33J investigated perfect fluid solutions using Bianchi type / space- 
time in scalar-tensor theory. Adhav et al. [33] obtained an exact solution of 
the vacuum Brans-Dicke field equations for the metric tensor of a spatially 
homogeneous and anisotropic model. Pradhan and Amirhashchi [35] inves- 
tigated anisotropic dark energy Bianchi type III model with variable EoS 
parameter in GR. Adhav et al. [3S] explored Bianchi type III cosmological 
model with negative constant deceleration parameter in Brans Dicke theory 
of gravity in the presence of perfect fluid. 

In this paper, we focuss our attention to explore the solutions of anisotropic 
dark energy Bianchi type III cosmological model in the context of BD theory 
of gravity. We find the solutions using the assumption of constant decelera- 
tion parameter and power law relation between <fi and a. The paper is orga- 
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nized as follows: A brief introduction of the field equations in BD theory of 
gravity is given in section 2. In section 3, the solutions of the field equations 
for Bianchi types III spacetime are found. Some physical and kinematical 
parameters are also evaluated for the solutions. Singularity analysis is given 
in section 4. A brief summary is given in the last section. 



2 Some Basics of Brans Dicke theory of Grav- 
ity 

The line element for the spatially homogeneous and anisotropic Bianchi type 
77/ spacetime is given by 

ds 2 = dt 2 - A 2 (t)dx 2 - e- 2mx B 2 (t)dy 2 - C 2 (t)dz 2 , (1) 

where A, B and C are cosmic scale factors and m is a positive constant. The 
energy momentum tensor for anisotropic dark energy is given by 

Tj = diag[p, -p x , -p y , -p z ] = diag[l, -u x , -u y , -u) x \p, (2) 

where p is the energy density of the fluid while p x , p y and p z are the pres- 
sures on the x, y and z axes respectively. Here u is EoS parameter of the 
fluid with no deviation and cu x , u y and ui z are the EoS parameters in the 
directions of x, y and z axes respectively. The energy momentum tensor can 
be parameterized as 

Tj = diag[l, -u, -(co + 7), -(u + S)]p. (3) 

For the sake of simplicity, we choose u x = 00 and the skewness parameters 7 
and S are the deviations from u on y and z axes respectively. 
The Brans-Dicke field equations are 

Rio - ^Rgij - ^{<t>,4,j - \gij4>,k4> ,k ) - ^(<P ; ij - giP4>) = (4) 



and 

ru — ^ k — 

3 + 2ro' 



^ = <^ = ^^> (5) 
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where w is a dimensionless coupling constant. For Bianchi type III space- 
time, the field equations take the form 

AB BC CA A 2 2 [ <j) ) fi B C> 0' [) 
B C BC 4> ©i 2 O C 87ro;p 
B + C + BC + + 2> } + ( } 

C i w 2 C A 87r(o; + 7)p 

C + I + O4 + + 2" ( } + ^C + A ] ~ 4> ' () 

A B m 2 AB I w <j> 2 <j> A B 8ti(uj + S)p 

a + b'a^ + ab + ^ + J { ^ + ^ { a + b ] - 4> ' {) 

Also, the 01-component can be written in the following form 

A B 

3 -5 = o. (io) 

Integrating this equation, we obtain 

B = Cl A, (11) 

where c\ is an integration constant. Without loss of any generality, we take 
C\ — 1. Using Eq. (JSJ), we get 

1 i,A B C\ 8tt(1 - 3a; - 8 - j)p 

* + «A + B + C )= ^(3 + 2^) (12) 

Now we define some physical parameters before solving the field equa- 
tions. 

The average scale factor a and the volume scale factor V are defined as 
a = \fAiC, V = a 3 = A 2 C. (13) 
The generalized mean Hubble parameter H is given in the form 

H = ^(H 1 + H 2 + H 3 ), (14) 

where Hi = ^ = H 2 , H 3 = ^ are the directional Hubble parameters in the 
directions of x, y and z axis respectively. Using Eqs.( TTB"|) and (113]) . we obtain 

H=\y / = \{H l + H 2 + H,) = -. (15) 
S V 6 a 
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The expansion scalar 9 and shear scalar a are defined as follows 

= »»=2^ + |, (16) 

° 2 - (17) 

where 

1 1 

<V = ^iu^-dK + u> v;a h%) - -9hp,, (18) 

— g^v ~ u^Uy is the projection tensor while u M = y^oo(l, 0, 0, 0) is the 
four-velocity in co-moving coordinates. The mean anisotropy parameter A m 
is defined as 

where AH t = H t - H, (i = 1, 2, 3). 



3 Solution of the Field Equations 

Subtracting Eq.flJJ) and Eq.(jHJ), we get 

B A C,B A s 6.B A x 8%jp . . 

Using Eq. flTTj) . this equation gives 7 = 0. Thus Eqs.fl&HH]) and Eq.( fl2|) reduce 
to 

,A. 9 2CA m 2 vj ,6. ~ 6 .A C\ 8np . . 

</ + M7-iy+h + c )= f' (21) 

(7 A CA 6 w f d> s0 6.C A. 8nup , . 

2A A 2 m 2 w0 2 2<M_ 87r(o; + (5)p 

•■ } .2i4 C\ 8tt(1-3o;-5)p 

* + «T + C ) = ^(3 + 2 CT ) (24 > 

Integration after subtracting Eq. (1221) and Eq. fT23l yields 

A C 1 f,m 2 8n5p^ 1 a2 „ , X 



where A is an integration constant. The integral term in this equation van- 
ishes for 

'-«& « 26 > 

Using Eq.()26]) in Eq.@5), it follows that 

A \ r dt 

-=c 2 e Xj ^, (27) 

where a 3 = A 2 C and C2 is an integration constant. Here we use the power 
law assumption to solve the integral part in the above equations. The power 
law relation between scale factor a and scalar field cp has already been used 
by Johri and Desikan [37] in the context of Robertson Walker Brans-Dicke 
models. Thus the power law relation between (p and a, i.e. (p °c a m , where m 
is any integer, implies that 

(f) = ba m , (28) 

where b is the constant of proportionality. The deceleration parameter q in 
cosmology is the measure of the cosmic acceleration of the universe expansion 
and is defined as 

<! = --■ 29 

cr 

It is mentioned here that q was supposed to be positive initially but recent ob- 
servations from the supernova experiments suggest that it is negative. Thus 
the behavior of the universe models depend upon the sign of q. The positive 
deceleration parameter corresponds to a decelerating model while the nega- 
tive value provides inflation. We also use a well-known relation between 
the average Hubble parameter H and average scale factor a given as 

H = la~ n , (30) 

where I > and n > 0. This is an important relation because it gives the 
constant value of the deceleration parameter. From Eq.f JT5|) and Eq.f l30|) . we 

get 

a = la- n+1 . (31) 

Using this value, we find that deceleration parameter is constant, i.e. q = 
n — 1. Integrating Eq. (j3ip . it follows that 

a = (nit + ki)», n^O (32) 
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and 

a = k 2 exp(Zt), n = 0, (33) 

where k\ and k 2 are constants of integration. Thus we obtain two values 
of the average scale factor that correspond to two different models of the 
universe. 

3.1 Dark Energy Model of the Universe when 

Now we discuss the model of universe when n^O, i.e., a = (nit + For 
this model, cf> becomes 

= b(nlt + hi)~%. (34) 

Using this value of 4> in Eq. (l2"T|) . the metric coefficients A, B and C turn out 
to be 



n-l 



A = B = c 2 Hnlt + k 1 )^eM X{n l^y i) n ), n*l (35) 



n-l 



3o(n — 1) 

The directional Hubble parameters Hi (i = 1,2, 3) take the form 

#i = #2 = ^ + n (37) 

nZi + fci 3b(nlt + h)i 

H 3 = ' (38) 

+ 3b(nlt + k 1 )^ 

The mean generalized Hubble parameter becomes 

# = , - , (39) 
n/t + fci 

while the volume scale factor turns out to be 

V=(nlt + k 1 )n. (40) 
The expansion scalar 9 and shear scalar a take the form 

9 = ^h: (41) 

a 2 = ^ T (42) 

36 2 (n/t + fc^n 
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The mean anisotropy parameter A m becomes 

2A 2 2 
A m = —(nlt + h) 2 —. (43) 

Moreover, Eqs.(j2TI)-(j23j) take the form 
8np , 9/ . . , , . 9 r A 2 / 2 m 2 — 2X(nlt + kA^ . . . , , . 2 

_^ = _/2 (3 + 2ro)(n ^ +A , i) -2_ [ __ + _^ exp( K _ V ^ nU+h y-^ 

(44) 

- ^ = / 2 [(3 + 2n + 2w){nlt + A4)- 2 + ^(n/t + h)^}, (45) 
3tr 

8n(uj + S)p , 9/ „ „ „ w , , , 9 r A 2 / 2 m 2 — 2X(nlt + /ci) 1 ^ 1 . , , , , , 2 

1 T = / 2 (3+2n+2w)Ht+A; 1 )- 2 +[ w r exp ^ )}(nlt+h ~. 

36 z C2 3 36(n — 1) 



(46) 

3.2 Dark Energy Model of the Universe when n = 0. 

The average scale factor for this model of the universe is a = k 2 exp(lt) and 
hence takes the form 

b 

= exp(-2Zt). (47) 
k 2 

Inserting this value of in Eq. (1271) . the metric coefficients A, B and C 
become 

A = B = c 2 h 2 exp(lt)eM- Xe ^~ lt) l (48) 

C = c 2 -h 2 eMlt)eM 2Xe ^ H) ]- (49) 
The directional Hubble parameters Hi become 

-3 - 

while the mean generalized Hubble parameter and the volume scale factor 
turn out to be 

H = l, V = k 2 3 exp(3lt). (52) 
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The expansion scalar 6 and shear scalar a take the form 



A 2 exp(-2/t) 



9 = 31, a 2 = (53) 

3b 2 k 2 2 V ; 

The mean anisotropy parameter A m for this model become 

= 2A 2 exp(-2/t) 

9b 2 l 2 k 2 2 V } 

For this exponential model of universe, Eqs.f l2ip -(l23l) take the form 

8np , 9/ , A 2 exp(— 2lt) m 2 , , 2A exp( — It) . , , 
^ = _ Z a( 3 + 2n7 p 2 ; ^-^exp -2Zt+ ^ >- ), 55 

gTT^p ,2/o, x , A 2 exp(-2/t) 
-^ = /(3 + 2ro)+ 3^ 2 2 ' (56) 
87r(a; + J)p 2 Q _ A 2 exp(-2/t) m 2 2Aexp(-/t) 
3b 2 k 2 c 2 zk 2 

(57) 



4 Singularity Analysis 

The Riemann tensor is a useful tool to determine whether a singularity is 
essential or coordinate. If the curvature becomes infinite at a certain point, 
then the singularity is essential. We can construct different scalars from the 
Riemann tensor and thus it can be verified whether they become infinite 
somewhere or not. Infinite many scalars can be constructed from the Rie- 
mann tensor, however, symmetry considerations can be used to show that 
there are only a finite number of independent scalars. All others can be ex- 
pressed in terms of these. In a four- dimensional Riemann spacetime, there 
are only 14 independent curvature invariants. Some of these are 

p p „ab p p p pab p p pabcd p pab pcd 

ti\ — ti — g JXab, til — tiabti , M-3 — ttabcd-ti , -"4 — K cd K ab . 

Here we give the analysis for the first invariant commonly known as the Ricci 
scalar for both models. 



10 



For model of universe with power law expansion, we can write Ricci scalar 
p(6 - 3n) A 2 / 2 m 2 

while for exponential model, it is given by 

A 2 m 2 
R = -2[6/ 2 + - 5 — ]■ (59) 

1 36 2 a 2 C2 § a 2 exp (2Aa ) J ^ ' 

Both of these models show that singularity occurs at a = 0. 



5 Concluding Remarks 

This paper is devoted to explore the solutions of Bianchi type III cosmo- 
logical models in Brans Dicke theory of gravitation in the background of 
anisotropic dark energy. We use the power law relation between (f> and a to 
find the solution. The assumption of constant deceleration parameter leads 
to two models of universe, i.e. power law model and exponential model. 
Some important cosmological physical parameters for the solutions such as 
expansion scalar 9, shear scalar a 2 , mean anisotropy parameter and average 
Hubble parameter are evaluated. 

First we discuss power law model of the universe. This model corresponds 
to n 7^ with average scale factor a = (nit + ki)~ . It has a point singularity 
at t = t s = — ^j. The physical parameters H±, H 2 , H 3 and H are all infinite 
at this point but the volume scale factor vanishes here. The metric functions 
A, B and C vanish at this point of singularity. Thus, it is concluded from 
these observations that the model starts its expansion with zero volume at 
t — t s and it continues to expand for < n < 1. 

The exponential model of the universe corresponds to n = with average 
scale factor a = k 2 exp(lt). It is non-singular because exponential function 
is never zero and hence there does not exist any physical singularity for this 
model. The physical parameters Hi, H 2 , H3 are all finite for all finite values 
of t. The mean generalized Hubble parameter H is constant while metric 
functions A, B and C do not vanish for this model. The volume scale factor 
increases exponentially with time which indicates that the universe starts its 
expansion with zero volume from infinite past. 
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The isotropy condition, i.e., — e > as t — > oo, is also satisfied in each 

case. It is mentioned here that the behavior of these physical parameters 
is consistent with the results already obtained in GR p]. The variable EoS 
parameter u for both models turn out to be 

Z 2 [(3 + 2n + 2w)(nlt + hi)' 2 + ^(nlt + h)'^) 

u — in ' 

/ 2 (3 + 2w){nlt + h)^ + [|£ + ^4 exp( - 2A g n + ^ " )){nlt + fc x )-S 

/ 2 (3 + 2^) + A2e 3 ^- 2 ^ 



/ 2 (3 + 2 W ) + + ^ exp(-2Zt + ) • 

Both of these equations suggest that at t = 0, oj has a positive value which 
indicates that the universe was matter dominated in its early phase of its 
existence. At t — > oo, the value of u turns out to be zero which indicate that 
the pressure of the universe vanishes at that epoch. 
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